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In the tight-binding theory, momentum matrix elements (MMEs) needed to calculate the optical 
properties are normally computed using a formulation based on the gradient of the Hamiltonian in 
the k space. We demonstrate the inadequacy of this formulation by considering the case of zigzag 
graphene nanoribbons. We show that one obtains wrong values of MMEs, in violation of the well- 
("^ ' known selection rules, if the unit cell chosen in the calculations does not incorporate the symmetries 

^vj , of the bulk. This is in spite of the fact that the band structure is insensitive to the choice of the 

unit cell. We substantiate our results based on group-theoretic arguments. Our observations will 
open an avenue for proper formulation of MMEs. 
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I. INTRODUCTION 



,JC^ . Tight-binding theory is one of the conceptually simplest, and widely used theories of the electronic structure of 
molecules, clusters, and solids. And indeed, the newly emerging field of graphene [I'j and related materials 0-ly| such 

QJ , as graphene nanoribbons (GNRs) (THIOJ has seen widespread use of the tight-binding approach for calculation of 

G ■ their electronic structure. |l ll4l5| and optical properties. [16] However, because of the unknown nature of the basis set 

• , associated with the tight-binding approach, the calculation of MMEs needed for computing the optical properties 

C^ ■ has always been a matter of debate. Blount |l7l| argued that the momentum operator for a periodic system can be 

G ' represented as 

^': p = ^VkTf (k), (1) 
o : 

CJ , where nio is the free electron mass, and Vki?(k), represents the gradient of the Hamiltonian in the k space. This 

'""'■ expression was used by Dresselhaus and Dresselhaus[i8| as well as Smith[l9| to perform early calculations of the 

_^ ' optical properties of solids using the tight-binding approach. Based upon generalized Hellmann-Feynman theorem, 

^ Lew Yan Voon and Ram- Mohan, ^2Qj argued that Eq. [l]is indeed the correct representation of the momentum operator 

\jQ for calculating optical properties. Adopting a gauge-invariant approach within the tight-binding formalism, Graf and 

^^ Vogel[21] obtained results in agreement with the work of Lew Yan Voon and Ram- Mohan, [20]. Cruz et aZ.[22| also 

^~^ indicated that, Eq. [T] leads to the correct computation of optical matrix elements. However, in a recent analysis, 

^~~l , Pedersen et aZ.|23| pointed out that Eq. [l]is incomplete, in that it does not contain the contribution of intra-atomic 

t~^ ' matrix elements. According to Pedersen et aZ.,^23| the MMEs (c(k)|p|w(k)) between the valence band states (|t;(k))) 

^^ ' and the conduction band states (|c(k))) is given by 

o ■ 

(c(k)|p|z;(k)) = ^^C:^(k)C„„(k)Vk(/3,k|iJ|a,k) 



X 

S' "^ n 



H' -K!I^^£i^^^^^^^C:,(k)a„(k)d,„, (2) 



a/3 



where, the valence band eigen state |f (k)) is expressed as 

|«(k))=^a„(k)|a,k); |a,k) = -^^e^''-^|a,R), (3) 



^ R 
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and a similar expression holds for the conduction band states |c(k)). In the equations above, R denotes a lattice 
vector, N is the total number of unit cells in the system, \a, R) is the a-th atomic orbital located in the unit cell at 
position R, -Bc,k(-E'«,k) is the energy eigen value of the conduction (valence) band and d^^ = (/3, 0|r|a, 0) is the matrix 
element of the position operator r with respect to the reference unit cell. Note that the second term on the right hand 
side (r.h.s.) of Eq. [21 called intra-atomic contribution[23], is the extra term as compared to Eq. [1] Sandu[24] further 
examined the issue and essentially agreed with the analysis of Pedersen et aZ.|23l| 

Recently, while developing a correlated electron approach for computing optical properties of GNRs[25l|, we discov- 
ered that the optical matrix elements of zigzag GNRs (ZGNRs) computed using Eqs. [1] or [5] were crucially depend 
on the choice of the unit cell. While the band structures of the ZGNRs in question, as expected, were found to be 
independent of the nature of the unit cell, however, correct values of optical matrix elements were not obtained unless 
the chosen unit cell also incorporated the point-group symmetry of the bulk. This result is counter-intuitive, and 
surprising, because normally we believe that the computed physical quantities for bulk systems should be independent 
of the choice of the unit cell. This issue is particularly important for the case of ZGNRs for which there is a certain 
ambiguity in the choice of unit cells, as compared to the case of armchair GNRs (AGNRs) for which a unique choice 
of the unit cell exists. To the best of our knowledge, this unit cell dependence of the formalism based on Eqs. [l]or 
[21 has not been reported earlier, therefore, here we aim to elaborate our findings, and to analyze our results, based 
upon group theoretic arguments. 

Remainder of this paper is organized as follows. In the next section we discuss the theoretical aspects of this work. 
In particular, we analyse the nature of unit cells chosen for the calculations from the point of view of their point 
group symmetries, and, based on their irreducible representations, deduce the optical selection rules. In section IIIII 
we present and analyse our results. In particular, we find that the results obtained are fully consistent with the 
optical selection rules deduced in section [TT] Furthermore, we also support our arguments by means of finite cluster 
calculations. Finally, in section [TVl we present our conclusions. 

II. THEORY 

In the present work, we consider a nearest-neighbor tight-binding (TB) model for GNRs, with zero site energies 

^= Y. t'Ui4aCja+c],C,,), (4) 

where (ij) implies nearest neighbors (NN) , c]^ creates an electron of spin a on the pz orbital of carbon atom i (assum- 
ing that the ribbon lies in the a;y-plane, with the a;-axis being the periodicity direction), and tij is the corresponding 
hopping matrix element. In order to obtain the band structure and the corresponding Bloch orbitals, the TB Hamil- 
tonian of Eq. [3] is Fourier transformed, and the corresponding matrix elements for the one dimensional (ID) system 
under consideration are obtained as 

H^Ak)^ E e'^'""*^(«a),,(0) (5) 

n— — oo 

where i{na) represents the i— th orbital of the unit cell located at position na, n being an integer and a is the lattice 
constant, j(0) represents the j— th orbital of the reference unit cell, and ii(na) j(o) is the corresponding hopping element 
which is non-zero only for the NN sites. The Hamiltonian obtained from Eq. [S] is diagonalized at different fc— points 
to obtain the band structure and the corresponding Bloch orbitals in the ID Brillouin zone. In order to compute the 
optical absorption spectrum within the TB model, the MMEs (c(k)|p|u(k)) need to be computed. For the purpose, 
we have used the formula proposed by Pedersen et a/.[23| (c/. Eq. [5]), as against the original approach of Blount [1^ 
embodied in Eq. [1] In order to ensure the correctness of the approaches, we performed calculations of the matrix 
element (c(k)|p2;|w(k))needed to compute the absorption spectrum for the light polarized along the x— direction, using 
both Eqs. [1] and[21 and found only quantitative differences. However, Eq. [21 is more general, and can also be used to 
compute the matrix element (c(k)|pj,|i;(k)) required for calculating the absorption spectrum for the y-polarized light, 
by setting the first term on its r.h.s. to zero, because for a ID system periodic along the x direction, the Hamiltonian 
has no ky dependence. 

For the present ID systems, we used the optical matrix elements to compute the optical absorption spectrum for 
the x-polarized (j/-polarized) photons in form of the corresponding components of the imaginary part of the dielectric 
constant tensor, i.e., exx {^yy i^^)) , using the standard formula 



v^c 



■K j a 



|(c(fc)bXA:))p 



^"^^'^^ /_./.{(i.™(fc)-M^+7^}i^.Vfc)''' ^'^ 




Figure 1: ZGNR-2, along with the two possible choices of unit cells which can generate it. The x-axis is taken along the 
periodicity direction, y— axis is perpendicular to it within the plane of the system, while the z-axis is perpendicular to that 
plane. Symmetry groups of cells (a) and (b) are C2v and d, respectively. Based upon the higher symmetry group of cell (a) 
we call it the symmetric cell, and (b) the asymmetric cell. The numbering of the atoms is consistent with the Hamiltonian 
matrices given in Eqs. [7]and|8l 



where i denotes the Cartesian direction in question, u represents the angular frequency of the incident radiation, 
Ecv{k) — Ec^k — Ey^k, 7 is the line width, while C includes rest of the constants. We have set C = 1 in all the cases 
to obtain the absorption spectra in arbitrary units. 

We will consider the optical matrix elements of various ZGNRs, characterized by their width parameter Nz, which 
is nothing but the number of zigzag lines across the width of the ribbon. In short, we will denote a ZGNR of width Nz 
as 7iGNK-Nz- We first consider the narrowest such ribbon ZGNR-2, shown in Fig. [1] which has four carbon atoms 
per unit cell and can be generated by periodically repeating either of the two different unit cells depicted in the same 
figure. Even though both types of unit cells lead to the same GNR in the bulk limit, the point group symmetries of 
the two unit cells are different. The unit cell of Fig. [T^ has C2v symmetry with the symmetry operators (in addition 
to the identity): (i) rotation by 180° about the a;-axis, (ii) reflection about the xy plane, and (iii) reflection about the 
xz plane, leading to four irreducible representations (irreps) Ai, A2, -Bi, and ,62. [ST] As per dipole selection rules of 
the C2V point group, z-polarized radiation will couple bands with Ai symmetry while the y— polarized radiation will 
couple Ai states to Bi states. [23| On the other hand, the unit cell of Fig. [Ud, has Ci symmetry with only the identity 
and the inversion being the symmetry operators. This group has irreps Ag and Au- and dipole selection rules allow 
x, y, and z, polarized (and mixtures thereof) radiation to cause optical transitions between Ag and Au bands. Thus, 
the symmetry analysis predicts that for the C2V unit-cell (Fig. [T^), the optical matrix elements will have either x- or 
y— non-zero components, while for the C,; case (Fig. ^p) both the components will be simultaneously nonzero. 

Using Eq. [31 and assuming that the NN hopping matrix element is t, the Hamiltonian matrix for the symmetric 
unit cell (c/. Fig. [T^) is obtained to be 



Hs{k) 



HAs(k) 



1 t(l + e-*'=°) 
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""ig- Eb) the corresponding matrix is 






/ i(l + e-*'''") 
i(l + e^^"-) t 
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t(l + e- 
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(7) 



(8) 



In order to compute the optical matrix elements (c/. Eq. [2), we need the matrix elements of the Vk^(k) operator, 
which for the present ID case, can be obtained easily by taking the derivative of each element of Hs{k)/HAs{k), with 
respect to k. Furthermore, the d matrix (c/. Eq. [5]) needed for the purpose is taken to be diagonal, with its elements 
being the Cartesian coordinates of various atoms of the reference unit cell. 

Next; we demonstrate that that the Hamiltonian matrices of Eqs. [7] and [S] are unitarily equivalent. Let us consider 
the fc-dependent unitary transformation 



U{k) 



10 

10 

1 

e 



(9) 



ika 



It is straightforward to show that Hs{k) 
will yield the identical band structure. 



U{k)^HAs{k)U{k). Therefore, it is obvious that both Hs{k) and HAs{k), 



Figure 2: Band structure of ZGNR-2 (cf. Fig. [ij obtained by diagonalizing either Hs{k) or HAs(k) {cf. Eqs. [7]and[8ll with 
t = —2.6 eV. Valence bands are labeled as vi/v2 and conduction bands as Ci/c2. 
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Table I: Irreps of the Bloch orbitals of ZGNR-2 at A: = and k 
unit cell. The Bloch orbitals are presented in the appendix. 
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III. RESULTS AND DISCUSSION 

Here we present and discuss our numerical results in two subsections, organized as follows: the first subsection 
contains the results for various ribbons in the bulk limit (infinite length), while in the second one the results are 
presented for finite-clusters of increasing sizes for the two choices of the unit cell. As shown below, the finite-cluster 
calculations help us in obtaining a better understanding of the results for the infinitely long ribbons. 



Infinite Ribbons 



The band structure of ZGNR-2 obtained by diagonalizing either Hs{k) or HAs{k) is identical, and is presented in 
Fig. [21 for t = —2.6 eV. The irreps of the corresponding Bloch orbitals at fc = and k = n/a, for the two choices of 
the unit cell are listed in Table U and can be easily deduced from the Bloch orbitals given in the appendix. 

In Table |TT] we present the values of optical matrix elements between different bands for both the symmetric and 
asymmetric unit cell cases, and for both types of unit cells the calculations were performed with Eq. [5] It is a well- 
known fact that in ZGNRs with even values of Nz, and for light polarized along the periodicity direction (a;-axis), 
HOMO-LUMO transition is forbidden, as are those between several other bands due to symmetry related selection 
rules. ^SJ Combining the knowledge of the irreps of various Bloch orbitals {cf. Table |l| at different fc-points, with the 
dipole selection rules of point groups C21, discussed in section |lll we note from table |ll] that, consistent with these 
selection rules, the x— component of the optical matrix elements between bands ci — vi (HOMO-LUMO) vanishes 
only if the symmetric unit cell is considered for the system. From the same table it is also obvious this transition is 
allowed for the y— polarized radiation when the symmetric cell is used, consistent again with the selection rules of the 
C2V point group. While, with the asymmetric unit cell, both the x and y components of the optical matrix element 
are found to be non-zero for the transition, a result in agreement with the Ci symmetry of the asymmetric cell. Thus, 
we obtain different results for the polarization characteristics of the radiation for the ci — vi transition: (a) with the 
symmetric cell it is strictly y polarized, while (b) with the asymmetric cell it has both x and y components. As far 
as the C2 — vi (and ci — U2) transitions are concerned both the unit cells predict it to be x polarized, however, the 
magnitude of the optical matrix element obtained with the asymmetric unit cell is much smaller as compared to that 
with the symmetric cell. 



Figure 3: Components of the imaginary part of the frequency dependent dielectric constant tensor of ZGNR-2, calculated using 



J2), 



,{2)/ 



t = —2.6 eV: (a) ei;tc (a;) for the symmetric unit cell (solid line), and asymmetric unit cell (broken line), and (b) eyy (u) 
identical for both the unit cells. A line width of 0.05 eV was used to plot the spectrum. 




In Fig. [3] we present the optical absorption spectrum of ZGNR-2, for both the x— and j/-polarized radiation, 
calculated for the two choices of the unit cell. From the band structure of ZGNR-2 (c/. Fig. [2]) it is obvious that the 
joint density of states (JDOS) of the conduction and valence bands have van Hove (vH) singularities at points k = 
and k — 7r/a, because of the parallel bands. This implies that, the optical absorption spectrum for the a;-polarized 
light, will have peaks corresponding to gaps between those bands at points fc = and k — n/a, for which the transition 

(2) 

is allowed by the selection rules. Therefore, we expect two peaks in e^a; (w) for ZGNR-2: (a) first peak around 2.5 eV 
corresponding to ci — V2/C2 — vi allowed transitions at fc = tt/o, and (b) the second peak close to 11 eV because of 

the transition between the same bands ed k — 0. For the y-polarized radiation, based upon the selection rules and 

(2) 
the vH singularities we expect three peaks in eyy (uj): (a) first one near 5 eV due to the C2 — V2 transition at fc = tt/o, 

(c) second one close to 8 eV due to the ci — vi transition at A; = 0, and (d) the final one near 13 eV caused by the 

C2 — V2 transition also at fc = 0. 

An inspection of Figl5)D reveals that the calculated values of eyy (w) exhibit precisely the three peaks (the last peak 

near 13 eV being barely visible) described above, both, for the symmetric, as well as for the asymmetric, unit cells. 

(2) 

The fact that the results for €yy{uj) are same with both choices of the unit cell is fairly obvious because the optical 
matrix elements needed to compute eyy (uj) depend only on the y coordinates of the sites (c/. Eq. [2]), which are 
identical for both the unit cells. This is also obvious from the y-components of the optical matrix elements listed for 
various fc-points in Table HIl 

For exx{i^), however, the situation is different. The calculated spectrum for the symmetric unit cell is fully consistent 
with the vH singularity based analysis of the JDOS, with the two peaks precisely at the predicted locations in Fig. 
[3^. However, in the same figure, the spectrum computed with the asymmetric unit cell agrees with this picture only 
for the first peak, while it contains two higher energy peaks located near 8 eV and 13 eV, in complete disagreement 
with the correct spectrum. From Table |ll] it is obvious that the optical matrix elements for both the symmetric and 
asymmetric unit cells are identical at the point fc — tt/o, therefore, even with the asymmetric unit cell we get the 
correct description of the first peak in exx{oj). However, the peaks around 8 eV and 13 eV in the asymmetric cell 
spectrum correspond to ci — vi and C2 — V2 transitions, respectively, at fc = 0, which have become allowed because of 

the incorrectly calculated values of the corresponding optical matrix elements (c/. Table lll|) . Thus, the peak around 

(2) (2) 

8 eV is now present both in exx{i^) and eyy (uj) computed with the asymmetric cell, leading to absorption with mixed 

polarization, a wrong result as discussed earlier. Furthermore, the peak close to 11 eV is missing from the asymmetric 

cell spectrum because for that case transitions ci — W2/C2 — vi have become wrongly disallowed at fc = (c/. Table lll| . 

Thus we conclude that with the asymmetric unit cell, for ZGNR-2 while one obtains the correct description of the 

lowest peak in the absorption spectrum, the predictions for the higher energy peaks in the spectrum are completely 

wrong. 

In order to ensure that the results presented here for the ZGNR-2 are universal, and not just valid for this particular 

ZGNR, we also performed similar calculations for the ZGNR-12. The symmetric and asymmetric unit cells considered 

for this ZGNR are presented in Fig. IH and again the symmetric cell has C21, symmetry while the asymmetric one 

possesses Ci symmetry. For both types of unit cells we obtain the identical band structure (not presented here), which 

is in agreement with that presented by other authors. [lg| Because the number of bands, and, therefore, the number 

of possible optical matrix elements is quite large for this case, in Fig. \5\ we directly present the optical absorption 



Table II: Optical matrix elements of ZGNR-2 between the bands identified in Fig. [2] at various points in the Brillouin zone, 
with k expressed in the units of vr/a. Column heading "Sym" corresponds to the values for the symmetric unit cell (Fig. [T^) 
and "Asym" denotes the results obtained for the asymmetric unit cell (Fig. [ijj). 
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(a) (b) 

Figure 4: Unit cells considered for the ZGNR-12: (a) symmetric unit cell, and (b) asymmetric unit cell. 



spectra for the a;-polarized photons (eia; (i^)) of this ribbon computed using both types of unit cells. 

From the figure, a trend similar to the case of ZGNR-2 is observed as far as the agreement between the spectra 
obtained using symmetric and asymmetric unit cells are concerned. The two sets of spectra agree at lower energy 
peaks, but diverge completely from each in the higher energy region. As compared to the symmetric cell spectrum, 
the one computed with the asymmetric cell in the high energy region has many missing peaks as well as several new 
peaks. The reasons behind this are precisely the same as for the case of ZGNR-2, that optical matrix elements are 
being wrongly computed with the asymmetric cell, thus allowing disallowed transitions and vice verse, leading to a 
spectrum in complete disagreement with the published results. |16l] In order to ascertain whether this behavior was an 
artifact of the tight-binding model, we included electron-electron interactions by employing the Pariser-Parr-People 
(PPP) model Hamiltonian[28] to perform similar comparisons. Although the results of our PPP model based work 
will be discussed in detail elsewhere. |25]but we again found this unit-cell anomaly while computing the optical matrix 




Figure 5: e\cJ;{ijj) of ZGNR-12, computed using t = —2.6 eV, for the symmetric (solid line) and the asymmetric (broken line) 
unit cells. A line width of 0.01 eV was used to plot the spectrum. 



Table III: HOMO-LUMO gap Eg (in eV), x- and y— components of dipole transition matrix elements {L\x\H), and {L\y\H) 
between the HOMO and LUMO orbitals, respectively (in A units), for finite clusters containing A'' primitive cells of the 
symmetric and asymmetric types. 
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Figure 6: Finite clusters containing A'^ repeat units of: (a) Symmetric cell, and (b) Asymmetric cell 
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B. Finite-Cluster Calculations 

ZGNR-2, with hydrogenated edges is nothing but the polymer polyacene, which has so far not been synthesized. [26'] 
However, several oligomers of polyacene (oligoacenes) such as naphthalene, anthracene, tetracene, pentacene, and hex- 
acene are known to exist, and their properties have been studied extensively both theoretically and ex per imentally. |26l| 
The symmetry group of oligoacenes is D2h, and the optical transitions in it are either x- or j/-polarized.[26] The HOMO- 
LUMO transition in oligoacenes is always y-polarized, just as in ZGNR-2, while the dominant a;-polarized transition 
is among the orbitals HOMO(LUMO) and LUMO+n(HOMO-n), where n{> 1) is an integer which depends on the 
oligoacene in question. [26] Noting that, the edge hydrogenation does not change the point group of the system, optical 
selection rules of ZGNR-2 and polyacene will be identical. Just as an oligoacene with an infinite number of repeat 
units is polyacene, an infinite number of repeat units of either the symmetric cell or the asymmetric cell will lead to 
ZGNR-2. Therefore, in what follows, we investigate the optical properties of clusters of increasing sizes, consisting of 
finite number of unit cells of the two types, to understand their evolution towards the bulk ZGNR-2. The shape of 
the finite clusters for the two cases is shown in Fig. [S] From the figure it is obvious that, because of the dangling 
bonds at the edges, the point groups of the clusters are C2v and Ci, for the symmetric and the asymmetric cases, 
respectively. Furthermore, without the dangling bonds, the symmetry group of these clusters will be D2h, just as in 
the case of oligoacenes. 

Table ITTT] co ntains our results on the HOMO-LUMO gap. Eg, HOMO-LUMO dipole transition matrix elements 
(_L|x|_ff)/A^[2^, and {L\y\H) for the increasing cluster sizes. The transition dipole moments were computed using the 
methodology adopted in our earlier work. [26] From the table it is obvious that: (a) For both types of clusters, the gap 
Eg is slowly closing with increasing N, (b) {L\y\H) is virtually the same for both types of clusters, and (c) {L\x\H)/N 
is exactly zero for all N for the symmetric case, but it has significant non-zero values for the asymmetric case, and 
demonstrates a slow decrease with increasing N. This confirms the results of the previous section that both types of 
unit cells yield the same band structure in the bulk limit, but different values of the x components of the transition 
dipole moments. 



IV. CONCLUSIONS 



Li this work, we investigated the influence of the nature of the primitive cells of ZGNRs on their electronic structure 
and optical properties, as computed within the tight-binding model. Based upon the Bloch orbital based calculations 
for the bulk systems, as well on finite clusters of increasing sizes, we conclude that the choice of the unit cell is 



of crucial importance when it conies to evaluation of optical matrix elements within /c-space gradient formulation, 
although, it is inconsequential as far as the band structure is concerned. In other words, Hamiltonians which are 
unitarily equivalent, lead to different values of optical transition matrix elements. This result is surprising because 
the theory of the optical properties of a system should not depend on the choice of the unit cell adopted to describe it. 
Our work demonstrate the inadequacy of this routinely used formulation for evaluation of optical matrix elements. We 
further substantiated our findings based upon the group theoretical analysis of optical matrix elements, and conclude 
that a unit cell which does not have the symmetries of the bulk system, leads to erroneous values of optical matrix 
elements violating symmetry based selection rules. Therefore, in order to describe the optical properties correctly 
within the fc-space gradient formulation, it is important to choose a unit-cell whose symmetry is consistent with that 
of bulk. 
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Table IV: Bloch orbital coefficients of ZGNR-2 at two k points for the symmetric cell. For each orbital, basis functions are 
numbered 1 to 4, from top to bottom. For each coefficient, the real part is followed by its imaginary part. 



fc = 


k = Tr/a 


V2 


Vl 


Cl 


C2 


V2 


"1 


Cl 


C2 


0.44, 0.00 


-0.56, 0.00 


0.56, 0.00 


-0.44, 0.00 


0.00, 0.00 


-0.71, 0.00 


0.71, 0.00 


0.00, 0.00 


0.56, 0.00 


-0.44, 0.00 
0.44, 0.00 


-0.44, 0.00 


0.56, 0.00 


0.00, 0.71 


0.00, 0.00 


0.00, 0.00 


0.00, 0.71 


0.56, 0.00 


-0.44, 0.00 


-0.56, 0.00 


0.00, 0.71 


0.00, 0.00 


0.00, 0.00 


0.00, -0.71 


0.44, 0.00 


0.56, 0.00 


0.56, 0.00 


0.44, 0.00 


0.00, 0.00 


0.71, 0.00 


0.71, 0.00 


0.00, 0.00 



Table V: Bloch orbital coefficients of ZGNR-2 at two k points for the asymmetric cell. For each orbital, basis functions are 
numbered 1 to 4, from top to bottom. For each coefficient, the real part is followed by its imaginary part. 



fc = 


k = Tr/a 


V2 


^1 


Cl 


C2 


V2 


VI 


Cl 


C2 


0.44, 0.00 


-0.56, 0.00 


0.56, 0.00 


-0.44, 0.00 


0.00, 0.00 


0.71, 0.00 


-0.71, 0.00 


0.00, 0.00 


0.56, 0.00 


-0.44, 0.00 


-0.44, 0.00 


0.56, 0.00 


0.00, 0.71 


0.00, 0.00 


0.00, 0.00 


0.00, 0.71 


0.56, 0.00 


0.44, 0.00 


-0.44, 0.00 


-0.56, 0.00 


0.00, 0.71 


0.00, 0.00 


0.00, 0.00 


0.00, -0.71 


0.44, 0.00 


0.56, 0.00 


0.56, 0.00 


0.44, 0.00 


0.00, 0.00 


0.71, 0.00 


0.71, 0.00 


0.00, 0.00 



Appendix A: Bloch Orbitals of ZGNR-2 for the symmetric and asymmetric unit cells 



The Bloch orbitals of ZGNR-2, obtained with the nearest neighbor hopping matrix element of -2.6 eV is presented 
in table IIVI for symmetric cell and in table |V] for asymmetric cell. 



